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Abstract The Australia Telescope 20GHz (AT20G) survey is a large area 
■ (27r sr), sensitive (40niJy), high frequency (20 GHz) survey of the southern 

sky. The survey was conducted in two parts: an initial fast scanning survey, 



Q ■ and a series of more accurate follow-up observatio ns. The follow-up survey 
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catalog has been presented bv lMurphv et al. In this paper we discuss 



> 

^ ■ the hardware setup and scanning survey strategy as well as the production of 

the scanning survey catalog. 
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1 Introduction 

The radio-source population above 5 GHz has not been well studied, in most 
part due to the lack of large-area surveys. Large aperture radio telescopes typ- 
ically have fields of view of a few arc-minutes at high frequencies, making it 
extremely time-consuming to carry out such surveys. The high-frequency prop- 
erties of radio sources are important to our understanding of the nature and 
evolution of both radio galaxies and the universe itself. Radio galaxy evolu- 
tion can be explored via AGN, whilst the nature and evolution of the universe 
can be explored via studies of the cosmic microwave background (CMB). High 
resolution radio surveys are crucial for the removal of foreground radio sources 
(Ide Zotti ct al.^ l99i) for CMB missions such as Planck (Taubcr . 2005i) . Fermi 
(jRitz et all . l2009l) and the South Pole Telescope (|Stark et all Il998l ). where 



foreground point sources contaminate the measured CMB anisotropy at angu- 
lar scales of less than 30 arc-minutes. 

A Wide-Band Analog Correlator (WBAC) was built as a prototy pe for the 
Taiwanese Array for Microwave Background Anisotropy (AMiBA iLo et al 



20011 ). This was not used in the production version of the AMiBA project, 
due to the complexity needed to remove 1/f noise in a large correlator, but 
it was realised that the WBAC could be used on the ATCA for a large area 
high frequency survey on a timescale that was not prohibitive. Since the time 
needed to survey a region of sky to a particular flux density depends both 
on the size of the primary beam and the integration time per pointing, the 
8 GHz bandwidth offered by the WBAC corresponds to a 32 fold increase in 
observing speed for a survey of the southern sky, as compared to the 256 MHz 
bandwidth available on the ATCA in 2004. 



2 AT20G Hardware 

The AT20G survey was conducted using three of the six antennas that make 
up the ATCA, along with a custom Wide-Band Analog Correlator (WBAC). 
In this section we describe the antenna and correlator design. 



2.1 ATCA telescopes 



To minimize the selection against extended sources in the survey, the most 
compact ATCA configuration was used. Antennas 2, 3, and 4 of the ATCA 
were situated on stations W102, W104, and W106, giving two 30.6m (2-3 
and 3-4) baselines and a si ngle 61.2m (2 - 4) bas eline. The antenna design and 
construction is outlined in Prater et al.l (Il992f). and the 20 G Hz upgrade to 



the Co mpact Array is discussed in lWong fc Melatoa ( 20021 ) and lMoorev et al 
(|2008l) 
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2.2 WBAC design 
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Fig. 1 Correlator architecture for a single baseline. The multipliers are located at the 
positions marked by and the path length differences correspond to the delay spacing of 
At. The complete WBAC consists of three such correlators. 



In 2002/2003 the WBAC was used in a pilot survey as a two element 
3.4 GHz bandwidth interferometer at 18 GHz. The pilot survey covered ap- 
pro ximately 12 deg^ of the sky between —59° > S > —71°, for further details 



see 



Ricci et al 



()2004D . The wide-band analogue correlator was adapted and 
expanded to produce a 3 element 8 GHz bandwidth, dual orthogonal polariza- 
tion interferometer with a central frequency of 20 GHz, which was used for the 
main AT20G survey. 



2.2.1 Assembly 

The WBAC is a lag correlator in which the baseline correlation functions are 
calculated firstly, and the cross power spectrum secondly. The operation of 
a lag correla t or, an d comparison against an FX architecture are discussed in 
Holler et al ] (l2007h . He we restrict our discussion to details specific to the 



WBAC. The lag correlator is constructed from an array of 16 multipliers on 
each of the three baselines, with precise transmission line delays providing 
the required delay difference for each lag. The architecture chosen is shown in 
Figure [1] and comprises two wide bandwidth power splitters to replicate the 
input signals, which are then arranged as two parallel counter-propagating sets 
of signals. The correlators were designed to process a 0-12 GHz input band to 
avoid having to perform a final down-conversion of the 4-12 GHz IF band to 
baseband. At the intersection of each corresponding pair of transmission lines 
is placed a multiplier. To facilitate assembly and simplify any required rework 
the multipliers are mounted on small second level microwave circuit boards 
which are inserted into cutouts milled in the correlator boards at the required 
positions. 

The separation of the lag channels. At, is determined by the displacement 
of the multiplier boards as shown in figure [TJ and is designed to be a half 
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wavelength at 12 GHz (ie Nyquist), which is a delay of 4.2 x 10"'^ nsec. The 
sampling in time delay is not perfectly uniform due to manufacturing and ma- 
terial effects, and may also have a small frequency dependence. This deviation 
from ideal Nyquist sampling is stable in time and can be accurately measured 
so the spectrometer's performance is not significantly affected. Figure [2] shows 
the transform phases of the WBAC as a function of channel number, for the 
primary calibrator 1921—293. The deviation from a linear fit to the phases is 
also shown, and is typically 0.2 x 10~^ nsec, which is less than 5%. The non 
ideal sampling is overcome by additional data processing steps as outlined in 
section [S] In a digital correlator the delay spacings are able to be adjusted so 
that the delay spacings are all equal and would result in a flat line at Onsec 
for the right hand plot of [5] 



5 10 15 5 10 15 

multiplier (c) multiplier (d) 

Fig. 2 Left: The transform phases for the WBAC with reference to an arbitrary time for 
the primary cahbrator 1921—293. The spacing between the correlator lag channels is of the 
order 4.2 X 10"'^ nsec. Right: The deviation from a linear fit to the transform phases. The 
largest error is 0.7 X 10~^ nsec with a typical value of 0.2 X lO"'^ nsec. 



2.2.2 Data Acquisition 

Timing and data acquisition of the correlator system is controlled from a Linux 
based computer, running a modified version of the standard Compact Array 
correlator control software. High accuracy timing signals are taken from the 
observatory's 5 MHz reference signal which is derived from a hydrogen maser. 
The data from the ATCA antennas are transported via an optical fiber to the 
control building where it is then converted to an electric signal by an optical 
receiver. The three antennas each have a phase switch in the signal path just 
after the optical receivers. The phase switching is designed to be mutually 
orthogonal and take the form of a 350 kHz square wave and two 700 kHz 
square waves. The signals are correlated and sent to an analogue to digital 
converter (ADC) where they are accumulated for an integration cycle before 
being sent to the control computer where the data are archived. 
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2.2.3 WBAC Performance 

Stability of the WBAC system was measured by observing the variance in the 
integrated spectrum as a function of time using a laboratory noise source and 
amplifier providing nominal input levels. This variance reduced in accordance 
with the radiometer equation, as the inverse square root of the integration 
time. This measurement was only taken over 10 minutes but was adequate to 
verify stability for the current transit instrument, which operates with inte- 
gration time less than 100 ms. The stability of the WBAC is important for the 
flux calibration of the survey observations as it allows for the use of a single 
bandpass and flux calibrator. 




Fig. 3 Schematic representation of the combined antenna and analog correlator system. 



2.3 Combined System 

Figure 13] shows a block diagram of the analogue correlation system. The system 
is comprised of three main parts: The IF transmission system that transports 
the down-converted received signals to the central site; the correlators; and 
the data acquisition and processing system. The ATCA analogue correlator 
lacks a variable delay compensator so it can be used only with transit inter- 
ferometers. The static delays are adjusted so that the array is directed towards 
the meridian with the antennas in a East- West linear array. 

2.3.1 Performance 

Once the WBAC was connected to the telescopes further characterisation and 
alignment of the system was performed. We adjusted the static delays to max- 
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imize the fringe amplitudes at the time a point source passed through the 
primary beam. 

The fringe pattern produced by each lag channel can be Fourier inverted 
to give the entire system frequency response for each correlator channel. From 
this data each of the lags can be calibrated. FigureHfLeft) shows the recorded 
fringe pattern for a single lag within the correlator as a calibrator source was 
observed. In figure |4l[ Right) the Fourier inverted system band-passes for all 16 
lags in the correlator are superimposed on a linear scale. The characteristic 
saddle shape is dominated by the gain curve of the 4-12 GHz amplifiers used 
in the IF transmission system. 




Time(x54ms) freq (GHz) 



Fig. 4 ATCA-(-WBAC system response for AT calibrator 1921-293. Left: Delay beam 
fringe pattern for a single lag channel on baseline 2 — 3. Right: A Fourier transform of a 
lag channel gives the system bandpass. All 16 channels are superimposed. 

The fundamental determinant of performance is the sensitivity achieved 
with the system. Observations with the 8 GHz, three baseline system were 
begun in September 2004. During these observations an average rms sensitivity 
of approximately 10 mJy, using 80 ms integration time, was achieved. 



3 Scanning strategy 

As it was not possible to insert the geometric delays necessary to track a 
point on the sky at 8 GHz bandwidth, the survey was intentionally restricted 
to meridian observations. By scanning the telescopes close to the maximum 
rate of ~15°min~-'^ along the meridian, the natural rotation of the earth al- 
lowed the observation of interleaving tracks across the sky. The position of 
the antennas was recorded every '^50ms. This gave a 0.8arcmin spacing be- 
tween observations along each track. The telescope slew speed and integration 
times were fine tuned so that the tracks would produce closed paths (over a 
24 hour period), with each subsequent path being one half-power beam width 
(2.3arcmin) from the previous scan. Figure [S] shows the interleaved tracks that 
are formed when data from multiple days are combined. The closed path allows 
for the distinction between tracks (interleaves), which is important to ensure 
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Fig. 5 An example of a scanning interleave pattern for —40° > 5 > —50°. A total of 7 
interleaves were required to achieve an interleave spacing equal to the primary beam FWHM 
of 2.3 arcmin. Data that were flagged from interleave (red) near 14'' 12™, is not a terminal 
problem as the same area of sky is covered once by observations scanning in the opposite 
direction. 

that any data which are missed or flagged due to bad weather or instrumental 
failure is easily repeated. The separation of the tracks by a beam width means 
that each point on the sky falls within the FWHM of the primary beam of the 
telescope. Sky locations on the interior of the scan region are covered twice 
each by this scanning strategy, once in a northward track and once in a (dif- 
ferent) southward track. Towards the boundaries of the declination band the 
coverage becomes less uniform and points are either covered only once, at the 
extreme edges, or twice with a single track. 

The survey mode observations for any given declination band of the sky 
were carried out contiguously, with all the observations for a particular region 
of the sky being completed within a single observing period of typically 10 
days. Observing in this way takes advantage of the stability of the correlators. 
Table [1] shows when the observations for each declination band were taken. 
To make the best use of the available time, any observations that needed to 
be redone were collected and generally carried out towards the end of the 
allocated observing time. For small amount of bad data that occurred towards 
the beginning or end of a particular interleave, the patching was done during 
the change over to the next day's observation. Due to our meridian scanning 
technique and the randomness of bad weather or equipment failure, it was 
not always possible to recover all the data which were lost or flagged. When 
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patching, preference was given to regions with missing data, and then to those 
with large amounts of noise or bad weather. 
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Table 1 The observation dates for the scanning observations for each of the declination 
bands. Overlapping dates are due to patching data being observed at the end of each year's 
observations. tThese regions include the south celestial pole and as such, required a separate 
program to produce the maps. 



4 Survey Calibration 

As mentioned in section !^. 2. 11 the WBAC has lag channel spacings that are not 
regular. In the absence of noise, the irregular lag channel spacing will not affect 
the recovery of the observed spectrum. With noise present, noise correlation be- 
tween lag channels will increase as the spacing diverges from Nyquist. To avoid 
correlated noise signals, the calibration and data reduction of the scanning in- 
terleaves was carried out in the lag domai n. Calibration of the lag re sponse 
can be done using monochromat ic signals (Harris fc Zmuidzinasl 2001 ). or in 



situ using astronomical sources (iHoller et al.l . l2007h . The AT20G survey ha d 



already adopted the approach independently described bv lHoUer et al. ( 2007 ). 



It is necessary to be able to deduce the location of a source in delay space 
from the signal distribution of the 16 lag channels of the correlator. If the 
calibration is correct over the three correlators then we will derive a common 
delay when a source is within the antenna beam. The source's position on 
the sky is derived from the delay, the antenna elevation, and the UT of the 
observation. Once the correlators are calibrated this process can be inverted: 
that is, given a record of 16 lags over 3 baselines, it is possible to estimate, for 
any delay, the amplitude of a source at a particular location which is consistent 
with the record. The key to this process is the tracking calibration which gives 
the signatures for all the delays within the correlator's range. A source on the 
equator observed with a 30m baseline will have 75 recorded samples in the 
time it takes for the peak correlation to advance from one multiplier to the 
next. 
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4.1 Overview 

The calibration scheme used for the scanning interleave data has five main 
components, two of which are instrumental and three of which are environ- 
mental. The instrumental calibration components are the round trip phase 
information, and the array alignment correction. These are measured and cor- 
rected for on timescales of less than a minute. The environmental calibration 
components use astronomical sources to detect variations in the correlator de- 
lay response and make the required adjustments. Figure [5] shows an overview 
of the calibration contributions. 



Round Trip Phase 
Delays 



Array Alignment 
Delays 



Primary Calibrator 
Beam-Shape and Flux 



Instrumental Calibration 
<1 IMinute 



Secondary Calibrator 
Delays 



Environmental Calibration 
~1 Day 



Tertiary Calibrator 
Delays 



Calibration 
of 

Raw Data 



Environmental Calibration 
~4 Hours 



Fig. 6 An overview of tlie various delay and flux calibration contributions. 



As there were no standard programs available that could cope with the 
specific setup of the WBAC, it was necessary to create a suite of purpose built 
programs. The programs were written in FORTRAN 71 ^ and implemented as 



tasks within the MIRIAD environment (|Sault et al.l . Il995f) . The software design 



and implementation is detailed in an upcoming ATNF memo. 



4.2 Timing Calibration 

The WBAC consists of 3 correlators, each correlating 16 lag channels for a 
single baseline. The output from each of the correlators' lag channels can be 
represented by a single baseline single channel interferometer as depicted in 
figure [T] The following is a simplified derivation of the received signal from each 
of the correlators, where noise terms, phase differences, and local oscillator 
mixing effects have been neglected for clarity. 

The signal received from each of the amplifiers can be written as: 

Va = A{l^) COSpTTi/t] 
Vb = Biv) COs[27ri^(t + T + Tast)] 



1 and later upgraded to FORTRAN 90 
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Fig. 7 Schematic of a single baseline, single channel east-west interferometer. The astro- 
nomical delay Tast is a function of meridian angle 8, and baseline b. Each correlator within 
the WBAC compute s C(Tast) for 16 di fferent delays (r) resulting in the function C{T,Tast)- 
Figure adapted from iThompsonl lll999l ) 






where: A{i/) and B{i') describe the voltage bandpass of the ampUfiers, the data 
transmission path from antenna to the correlator, and the gain equalisation 
units at the correlator, r is the compensating delay inserted by the correlator, 
and Tast is the astronomical delay: 

Tast = - sin(6') 
c 

for an East- West baseline of length b observing a source at meridian angle 9. 
The correlator receives the two signals and multiplies them to form: 

VaVb = M{u)A{iy)B{v) cos[2myt] cos[27riy(t + r + Tast)] 

= M{l^)A{l^)B{iy)^{cOs[2TTI^{2t + T + Tast)] + COs[27rt/(r + Tast)]} 

where M{i/) is the bandpass of the multiplier and the result been sep- 
arated into a high and low frequency component. A low pass filter is used 
to remove the high frequency component with frequency Anvt, leaving only 
the low frequency component which is then integrated to give the correlation 
function: 

C{T,Tast)= / M{l^)A{l^)B{l^) C0s[2Tn^ {T + Tast)]diy 

where Vc is the low pass filter's central frequency and Av is its bandwidth. 
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Fig. 8 The correlator response to the caUbrator 1921—293 on baseUne 2—3. All 16 lag 
channels are shown and the channel number is indicated in the upper right of each plot. The 
abscissa is time since transit in units of 54ms. The ordinate is un-calibrated flux in units of 

Jy- 



The basic observable is the set {C{T,Tast)}- Primary calibrator observa- 
tions are done in a tracking mode where by the calibrator is tracked by the 
telescopes as it transits the meridian. During this time the set C{T,Tast) is 
recorded for range of delays Tast that is large enough to encompass the synthe- 
sised beam response as is demonstrated in Figure[5] Having the antennas track 
the calibrator allows the primary beam response of the antennas to be disen- 
tangled from the synthesised beam response of the correlator. Once C{T,Tast) 
has been measured, A{v)B(i')M{v) can be derived via Fourier inversion, and 
the correlator response can be treated as given in the above equation. The sep- 
aration between lag channels {At) is A/2 where A = 24mm and corresponds 
to 85 arcsec in RA for a 30m baseline. The separation in Tast is the integration 
period of 54 ms and represents 1 arcsec in RA. 

In addition to the primary calibrator observation, timing adjustments are 
required for environmental, and instrumental effects that vary on time scales 
of minutes to hours. The main environmental factor is the phase difference 
between the signals from two telescopes due to their slightly different trans- 
mission paths. The ATCA is equipped with a round trip phase tracking system 
that measures the delays introduced into the telescope to correlator transmis- 
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sion path due to thermal effects within the conduits containing the optic fibers. 
The round trip phase is sensitive to diurnal changes in temperature within the 
conduits as well as the airconditioning cycles within the telescopes and corre- 
lator room. The main instrumental factor is the array alignment. The ATCA 
array deviates slightly from true East- West, requiring an elevation dependent 
correction to the timing solution. 

Secondary and tertiary calibrators are used to refine the timing corrections 
that are not otherwise captured by the ATCA monitoring system. The number 
of calibrators within the ATCA calibrator catalog is high enough that, on 
average, 3-4 daily calibrators will fall within the primary beam of the array 
during the normal survey scanning observations. Since each part of the sky is 
covered by two separate interleaves each of these secondary calibrators appears 
in at least two different interleaves and thus provides a way to calibrate the 
inter-scan timing solutions. To increase the number of calibrators available, 
strong {S20 > 100 mJy, SNR > 5) sources found within the calibrated scanning 
observations were used to self calibrate the timing, by requiring phase closure 
between the three baselines. The number of available calibrators was increased 
from around 500 to nearly 3000 over the entire sky using this bootstrapping 
process. 

Once the follow-up observations of the AT20G survey had been completed 
and a list of confirmed sources had been created, a final round of calibration 
was done using the improved positions and fluxes of all the sources with SNR> 
5. It is from this final round of calibration that the catalog described in section 
Owas created. 



4.3 Flux Calibration 

In the calibration scheme outlined in the previous sections, the primary fo- 
cus is on timing calibration. Whilst the calibration tables themselves contain 
both delay and amplitude calibration measures, the stability of the correla- 
tors meant that the primary calibrator was enough to ensure that flux densities 
were as correct as could be expected from the relatively short observation time 
of each source. For a typical source in the interior regions of a map that has no 
flagged data, there will be 6 observations of 54 ms each — three in a northward 
scanning interleave and three in a southward scanning interleave. This gives a 
total integration time of 324 ms and a system equivalent flux density (SEFD) 
of 300 — 400 Jy, corresponding to an rms noise level of 10 — 14mJy/beam. 

In the time between the scanning survey observations and the follow-up 
confirmation observations the absolute flux density scale for a particular re- 
gions was not accurately known. However, the stability of the WBAC meant 
that the relative flux density scale within a declination band was constant. 
When scheduling the follow-up observations, the strongest sources were given 
higher preference so that, in the absence of clustering, they were observed be- 
fore the weaker sources. Once the follow-up observations were complete and 
the data had been reduced it was then possible to scale the scanning survey 
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flux densities to that of the follow-up observations. This re-scaling does not 
add any extra information to any of the sources within the follow-up cata- 
log, but it is useful when computing the completeness of the two surveys as 
discussed in section [TUl 
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Fig. 9 Ratio of follow-up survey catalog and scanning survey fluxes. The horizontal black 
lines represent the mean flux ratio for each declination band, the values of which are printed 
below the data. The red lines represent a la deviation from the mean. The mean flux ratios 
were used to bring the survey map flux density scale in line with the more reliable follow-up 
observations. The year of scanning observations is given at the top of each declination band. 



Figure [5] shows a comparison between the follow-up survey catalog (FSC) 
flux densities and the scanning survey flux densities as a function of decli- 
nation. The scatter in the flux ratio is due to the fact that the follow-up 
observations were much more sensitive and have better (u, v) coverage, allow- 
ing for a more accurate measure of the flux densities. The variation between 
regions that were observed contiguously is much less than that of regions that 
were observed in different years. The flux of the band-pass calibrator sets the 
flux density scale for each particular region. The flux of 1921—293 varies slowly 
enough that it possible to interpolate the observations of the AT calibrator 
monitoring project to obtain an accurate flux. 



5 Motivation for a Scanning Survey Catalog 

Initially the purpose of the scanning observations was to detect individual 
sources so that more accurate, follow-up observations could be made of the 
stronger sources. The amount of time between scanning observations and 
follow-up observations was typically only 2-3 weeks in order to minimize the 
effects of source variability. This required that the calibration and production 
of candidate source lists be done as quickly as possible. The lists of candidate 
sources obtained from the scanning data were used only as an intermediary 
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between scanning and follow-up observations. The fluxes and positions had 
low accuracy and early on in the survey there were many spurious candidate 
sources which decreased the efficiency of the follow-up observations. If the 
accuracy of the positions and fluxes of the candidate source lists could be 
improved, and the number of false positive sources reduced or at least well 
characterised, then the candidate source lists become a useful catalog in their 
own right. Once the AT20G follow-up observations were completed the cali- 
bration and map making process was able to be finalized and applied to each 
of the declination bands to produce more reliable candidate source lists. It also 
became possible to use the final AT20G (follow-up) survey catalog to gauge 
the effectiveness of different calibration and map making procedures. 

The remainder of this section is devoted to creating the best possible cat- 
alog of sources using the scanning survey data from the AT20G. To avoid 
confusio n, the sources t hat co nstitute the final release version of the AT20G 
catalog ( Murphv et al. . 20101 ) will be referred to as follow-up sources, or as 
belonging to the follow-up survey catalog (FSC). The sources that were iden- 
tified within the survey maps will be referred to as scanning survey sources or 
as coming from the scanning survey catalog (SSC). Although the FSC should 
logically be a sub-set of the SSC, this is not strictly true. Of the sources within 
the FSC 89% of them are also within the SSC. The 11% of sources are missing 
from SSC due to the flux and reliability cut off that was imposed - most of the 
missing sources have a signal to noise less than 5, while the stronger missing 
sources have a reliability which is < 50%. Further details and discussion given 
in section ISl 

A catalog produced from the scanning observations is able to cover a larger 
area of sky than the subsequent follow-up observations and to a lower limiting 
flux density. The reliability of the claimed sources within the scanning survey 
catalog will not be as high as if follow-up observations had been done, as there 
is no way to directly confirm a source's existence. A catalog which is more 
complete and yet less reliable than the follow-up survey catalog will be of 
particular use for the construction of foreground masks for CMB anisotropy 
missions. In the measurement of the CMB anisotropics, excluding a real source 
from the mask is far more detrimental than the inclusion of a noise spike, and 
thus the most complete catalog is the most useful, even when not all of the 
claimed sources are real. A secondary use of the scanning survey catalog is 
to identify sources that were excluded from the FSC due to observing time 
constraints or poor weather. If the reliability and completeness of the scanning 
survey catalog is characterised well enough, it is then possibl e to es timate the 
completeness of the FSC. As was reported by iMurphv et aP (|20inh . the FSC 
is 91% complete at 100 mJy and 79% complete at 50mJy. 



6 Imaging the sky 



Once a suitable calibration scheme had been developed as described in section 
SI it was then possible to combine the scanning interleaves to produce a map of 
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the sky. Due to the non-uniform spacing of the delays within the correlator it 
was not possible to create a map using standard Fourier inversion techniques 
as an accurate inverse transform is no longer guaranteed, and correlated noise 
signals will result. 

As each part of the sky falls within the primary beam of more than one 
observation (see figure [S]) it was necessary to create an auxiliary map of weights 
that would allow the maps to be normalised once all of the observations had 
been processed. Pixels were weighted according to their distance from the 
beam center and in proportion to the square of the primary power gain, and 
the weights were scaled so that a pixel that was observed at the beam center 
had a value of 1. By scaling the weights in this manner it was possible to 
use this auxiliary map to determine how well each part of the sky had been 
covered by the scanning observations. This sky coverage map was used to give 
quantitative confidence levels in the source detection stage (section [S]) , as well 
as serving as a way to determine the survey sky coverage (section [T]) . 

The SEFD, sky coverage and calibration effectiveness are all critical to 
producing maps with high signal to noise. Figure fTOl shows the rms noise level 
for each of the maps within the survey region. The typical rms noise is around 
lOmJy/beam, which is in agreement with that calculated in section [4.31 for a 
SEFD of 300 Jy. The region —30° > <5 > —50° has a higher rms noise value 
of around 14mJy which is consistent with the measured SEFD of 400 Jy. The 
variation in rms noise level seen in the region —60° > 6 > —70° around 
a = is'' is due to patching observations that were done using the incorrect 
interleave number and thus areas with poor data were not recovered. 

The resulting (flux density) map of the sky remains convolved with the 
telescope response pattern. A CLEAN algorithm was not able to be used to 
deconvolve the map as a Fourier inversion would not produce a map consistent 
with our data. A specialised source detection program was created that would 
match sources to the telescope's beam pattern as measured by the primary 
calibrator. Further details are discussed in section |8l 



7 Sky Coverage 

As mentioned in section[B]the procedure that created the sky maps also created 
a map of the pixel weights which, when correctly scaled, become a map of the 
sky coverage. Figure [TT] shows a region of sky and the corresponding coverage 
map. The map of the sky coverage reveals a lot about how the survey was 
conducted and what sorts of problems were encountered. The interior parts of 
the map are scanned at least twice with a separation between scans being of 
the order of days, and the scans are in opposite directions. Figure [TT] shows that 
the typical sky coverage is around 3 ±0.5 scans/pixel, with variation due to the 
way that the different scans overlap at each point on the sky. As the telescopes 
require some time to turn around at the northern and southern boundaries of 
each declination band, the spacing between observations decreases and the 
sky coverage increases dramatically toward the boundaries. The coverage at 
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Fig. 10 The corrected rms noise level for the survey maps. The rms noise level is consistent 
with an SEFD of 300-400 Jy. Variations within a declination band are due to variations in 
the survey coverage. 



the turn-around points of the scans can be as much as 20 scans/pixel. When 
data are flagged for poor quahty on the interior of the map, and no patching 
has been done to recover the data, there is a reduction in the sky coverage at 
that point. The boundary regions of the maps, whilst still having a twofold 
redundancy in coverage, have a scan separation of only a few seconds in time, 
and thus data flagged for a contiguous block of time removes scans in both 
directions. The small gap in the map at 14''44™20'* - 14°57'23" is a resuh 
of data that were flagged out due to the presence of strong radio emission 
from Jupiter. When a significant amount of data are lost due to bad weather 
or otherwise, an attempt is made to re-observe the data and the patching 
data generally overlaps the flagged data by a minute or more. The result in 
this over-patching is better coverage for some parts of the sky. This is also 
evident in figure [TT] where the coverage increases to almost 6 scans/pixel in 
some places. 

For the most part, variations in the sky coverage do not affect the quality 
of the scanning survey catalog as they are only small. However there were 
problems with the zenith (-30°18'46" for the ATCA), that caused a lot of 
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Fig. 11 Top: Convolved sky map for a small region of the sky. There is a small amount 
of sky missing at the top of the map due to flagged data. Three sources with flux densities 
of 222 mjy (left), 286 mjy (middle), and 284 mjy (right) are easily visible. Bottom: Sky 
coverage map for the same region. Areas in white have no data, whilst black areas have a 
coverage lower than 2 scans/pixel. The effect of flagged data on the northern boundary can 
be seen. The scan in green is an instance of patching data that overlaps with good data and 
thus increases the sky coverage for a small part of the sky. 

data to be flagged at the hardware level and were not able to be recovered. 
The large amount of flagged data around —32° declination causes the sky 
coverage to drop below 1 scan/pixel. The amount of sky that is affected is 
small but the areas of low coverage produce many spurious sources that make 
the detection of real sources in this area quite difficult and thus decreases the 
completeness of the survey. 

8 Source Detection 

The way in which candidate sources were detected using the scan data changed 
over the course of the AT20G observations. In 2004 the scan data were pro- 
cessed as time ordered data using a stream-based approach. The final candi- 
date source detection of 2007 was based on the scan maps discussed in the 
previous section. Intermediate years used a combined approach with the de- 
tection rate of the followup observations increasing year by year. In each case 
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it was necessary to have a fast turn around for the candidate source hsts, as 
fohow-up observations were typicahy scheduled within 2-3 weeks of scanning 
observations to reduce the effects of variabihty. The two detection methods 
are discussed below. 

8.1 Stream-Based Detection 

In the early years of the survey, before the intra-scan calibration methods had 
been refined, candidate source lists were required quickly, so that follow-up 
observations could be done. The following equation: 

^(Ta.t) =5]C(r,r„,t)Cob.(r) (1) 

where Cobs are the measured correlator data, and C{T,Tast) system response 
to a calibrator (Figure [S]) , was used to identify candidate sources whenever 
a signal to noise of 5 or more was reached. This method was relatively fast 
and easy to compute but was less sensitive to sources that were towards the 
half power point of the primary beam. Also, due to the lack of intra-scan 
calibration, sources that appeared in two separate data streams often had 
position and flux discrepancies. The problem of side-lobe and multiple de- 
tections was significant and many side-lobe or duplicate sources entered into 
the follow-up observations. Once a follow-up observation had been made such 
occurrences were easily detected. In 2004 the confirmation rate of the follow- 
up survey was as low as 50%, meaning that up to half of the observing time 
allocated to follow-up observations resulted in non-detections. In 2005 the 
detection scheme, calibration, and duplicate/side lobe rejection methods were 
improved. The problem of reduced sensitivity to sources that were at positions 
not scanned close to the beam center remained. During 2005 the calibration 
scheme was improved to such an extent that a map based detection method 
was comparable in reliability to the stream based method, and a correlation 
of the two results was used as a candidate source list. From 2006 onwards the 
map based detection method was reliable enough that it was used alone. 

8.2 Map-Based Detection 

As the final maps were still convolved with the antenna response, it was neces- 
sary to create a specialised source detection program (detailed in an upcoming 
ATNF memo). The program works in a very CLEAN-like manner: searching 
for the largest peaks in the map, fitting to a predetermined model, and remov- 
ing it from the map. The primary calibrator was used to determine the an- 
tenna response which makes up the template against which candidate sources 
are matched. The antenna response to the primary calibrator PKS 1921— 293 
is shown in figure 1121 Starting at the highest peak in the map an iterative 
process was run where by the template was matched to the peak, a fiux and 
position were recorded along with: a measure of how well the source response 
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matched the cahbrator template, the number of pixels that were within the 
fitting area, the number of blanked pixels within the same area, and the sky 
coverage at the source position. A region of the map around the peak was 
then blanked and the process continued until some minimum flux density was 
reached. The maps were produced in blocks of l''6™ in RA and between 10° 
and 15° in DEC. The extra six minutes of overlap in RA between the map 
and its neighbours avoids problems involved when detecting a source on the 
boundary of a map. 




Right Ascension (J3000) 



Fig. 12 The convolved image of primary calibrator 1921—293, which was used to create 
the template used in the detection of sources within the convolved sky maps. The lowest 
contours are (red/blue) ±0.25 Jy/beam with consecutive contours doubling up to 8 Jy/beam. 
Primary beam side-lobes can be seen outside of the contours at a few percent of the source 
flux of 11 Jy 

The searching algorithm does not take account of possible extended sources. 
For the 30 m survey spacing there is a 50% reduction in amplitude for a Gaus- 
sian source with size 45 arcsec, and negligible response for sources much larger 
than 1.5arcmin. Any sources that are extended will have underestimated flux 
density in these survey maps. Fornax A, a well known extended radio source 
at low frequencies, is present within the maps at a level of 50mJy. All of the 
extended emission is completely resolved out. A point source that falls within 
a beam width of a stronger source will not be detected in the source detection 
program as there is no allowance given for double or multiple sources, however 
such sources will become evident during follow-up observations. The follow-up 
observations have an increased (u, v) coverage and integration time, resulting 
in a better position and flux density measurement. The FSC does not suffer 
as badly from the problem of missing flux for extended sources, although a 
source with no bright core or hot spots will likely not be detected in the survey 
maps and thus not make it into the follow-up observing schedule and thus not 
be present in the FSC. 
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The list of candidate sources found in the maps was fihered to remove 
side-lobes of strong sources and sources that were detected twice due to their 
position at the edge of two different maps. The Galactic plane is a plentiful 
source of extended and close multiple objects of sufficient strength to confuse 
the source detection algorithm and are not well characterised in the maps. 
Candidate sources with a Galactic latitude of \b\ < 1.5° were thus removed 
from the list. Areas of the sky which are close to the declination boundaries, or 
which have flagged or missing data in one or more interleaves, will have a much 
higher noise than surrounding areas. Since the source detection algorithm uses 
only a single noise value for an entire Ihour by 10° — 15° map these regions 
will increase the overall noise level by only a small amount, but contribute 
many false sources with a flux more than five times the map rms noise (see 
Figure [T4t . Spurious sources were identified and removed via a filtering metric 
which is now discussed. 



8.3 Source Filtering 

In order to increase the reliability of the scanning survey catalog, a filter 
was used to remove sources that were likely to be spurious.. Lower frequency 
catalogs were used to guide the creation of the filter but were not in the 
actual source selection. Once the initial candidate list had been prepared it 
was matched against lower frequency sky surveys to determine what limits 
should be placed on the filtering process to minimise both the number of 
spurious sources included and real sources excluded. A combination of the 
Sydn ey University Molonglo Sky Survey (SUMSS version 2.1 of March 11 
2008. iMauch et all (|2003l )). and the Second epoc h Molonglo Gi il actic Plane 



Survey (MGPS-2 version 1.0 of August 15 2007, iMurphv et all (|2007t )) was 



used to identify sources within the candidate source lists that were likely to be 
real sou rces. In regions above a declination of -30°the NRAO VLA Sky Survey 
(NVSS, ICondon et aP (Il998h 'l catalog was used. Whilst it was not expected 



that the lower frequency catalog would be completely recove red, it was known 
from previous observations (in particular Ricci et al. ( 2004 )). that almost all 



sources detected at 20 GHz were previously known in the more sensitive lower 
frequency surveys. 

A source was considered to be identified if had a low frequency counterpart 
within 1 arcmin. Some follow-up survey catalog (FSC) sources were initially 
found to have no counterpart in the relevant catalog, but closer inspection of 
the survey maps revealed extended sources that were either not listed in the 
scanning survey catalog (the MGPS-2 catalog contains only the point sources), 
or that have a hot-spot that is some distance from the FSC position. Of the 
approximately 6000 s ources within the AT20G FSC only 2 5 were unidentifie d 
(jMurphv et al.l . l2010l ). This result is consistent with that of lRicci et"all ()2004h . 



who found low frequency counterparts to all 221 sources detected at 18GHz. 

The matching rate between the scanning survey catalog and the low fre- 
quency catalogs is shown in figure 1131 where a strong correlation between flux 
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Fig. 13 Low frequency matching rate as a function of corrected survey flux at 20 GHz for 
the unfiltered scanning survey catalog. The blue line shows the number of sources per flux 
bin, whilst the red line shows the fraction of sources matched per flux bin. The dotted line 
represents the random matching probability of 5%. 



density and matching rate can be seen. This correlation is expected since the 
number of spurious sources in the scanning survey catalog increases as the 
survey flux decreases. The anti-correlation with source count is expected for 
the same reason. 

Figure [T3] shows all of the sources that were detected in the survey maps 
with an SNR > 5. The raw scanning survey catalog shows features that can 
clearly be attributed to noise, such as the over-density of sources towards 
the boundary of each declination zone. Variations in the survey coverage 
are also a cause of spurious sources and can be seen around 18 hours in the 
-50° > 6 > -60° zone. The zone -30° > S > -40° contains the zenith as 
seen by the ATCA and thus has a region of poor coverage interior to the map. 
The overlapping of sources towards the zone boundaries is seen at a much 
lower level in the matched source plot, which can be attributed to sources 
that were detected in the overlapping regions of maps from different declina- 
tion zones. It can be seen from Figure [TJ] that matching with low frequency 
catalogs is an effective way to distinguish between real and spurious sources. 
The correlation with lower frequency surveys was used as a guide to deter- 
mine what combination of the quality control quantities would produce the 
largest fraction of confirmed detections. The low frequency surveys were used 
to determine what a real source would look like within the map; a source's 
presence in a lower frequency survey was not part of the filter, so as not to 
bias the final catalog. At sufficiently low 20 GHz flux densities the correla- 
tion between low frequency detection and foUowup confirmation is expected 
to break down. For a source to be detected at 50 mJy in the AT20G a spectral 
index of °4 g § > -1-1.3 would be required to avoid detection in the NVSS, 
SUMSS, or M GPS-2 surveys. This ra r e class of ultra inverte d spectrum sources 
is discussed in Massardi et al. ( 2008 ): Murphv et al. ( 2010h . 

With a matching radius of 60 arcsec, some matches will be du e to chance 



and no t indicative of any real identification between the two sources. iCondon et al 



(|l998f ) give a probability of 5% for finding an unrelated radio source within 
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60 arcsec of an NVSS source. The SUMSS and MGPS-2 catalogs have source 
densities of 30 sources per square degree, whilst the NVSS has 50 sources per 
square degree. Thus we expect that around 5% of the identifications within the 
middle plot of figure [T?] are spurious. Figure [T3] shows a minimum matching 
rate that is not zero and can be attributed to this effect. 
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Fig. 14 Top: All the sources that are found in the survey maps with a signal to noise 
greater than 5. Middle: Sources that are matched with a low-frequency source. Bottom: 
Sources that are not matched with a low-frequency source. The matching with the NVSS, 
SUMSS, and MGPS-2 catalogs is a good indicator of sources that are likely to be real. 



To distinguish between real sources and those that are likely to be due to 
poor coverage or high noise peaks, a filtering metric was established with the 
guidance of a genetic algorithm. It was found that the metrics that produced 
the greatest separation between the matched and unmatched source popula- 
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tions were variations on equation [5] Here 6*20 is the flux in Jy as measured in 
the survey map, is the cross correlation of the source and the cahbrator 
template, cov is the coverage as measured in the auxiliary map, and pix is the 
ratio of good to bad pixels that were within the blanking region around the 
source. 

^ = logio ('5'20-X^-Cov.pix) (2) 



Unmatched at low frequency Matched at low frequency 




Survey Flux (Jy) Survey Flux (Jy) 

Fig. 15 The filtering metric as a function of survey flux for the scanning survey catalog. 
Blue sources (right) are those that have a low frequency counterpart within a 1 arc-min 
radius, red sources (left) are those that do not. The horizontal line represents the point at 
which equal numbers of sources are matched and unmatched. 



Figure [T5] shows as a function of survey flux for sources within the scan- 
ning survey catalog, identified via matching with the low frequency catalogs. 
The horizontal line at -2 represents the point above which there are an equal 
number of matched and unmatched sources. A catalog constructed with only 
sources that have J- > —2 will have an overall matching rate of 50%. It is 
expected that nearly all the real sources at 20 GHz will have a low frequency 
counterpart, however if that counterpart is complex or extended then a more 
careful match between the catalogs will be required to identify counterparts. 
For this reason it is expected that some of the 20 GHz sources with J- > —2 
will remain unmatched despite the fact that they are real sources. Similarly 
not all the sources that are matched at low frequency are expected to be real. 
In figure [Wl riohi) there are 1758 sources below J' = —2, which is 4.7% of 
the 35598 unmatched sources in the same area of the plot. This is in agree- 
ment with the 5% random matching rate cited for the low frequency surveys 
previously. The matched and unmatched sources occupy the difi'erent, yet over- 
lapping, areas of {J^,S) space. Using the low frequency identification as a proxy 
for the detectability of a survey source, it is apparent that eq [2] is efficient at 
separating real sources from noise peaks. Figure [T51 is the same as ll5l with the 
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Fig. 16 The filtering metric as a function of survey flux for the scanning survey catalog. 
Blue sources (right) are those that have a follow-up survey catalog counterpart within a 1 
arc-min radius, red sources (left) are those that do not. The horizontal line represents the 
point at which equal numbers of sources are matched and unmatched at low frequency. 



sources matched with the FSC. The same separation of sources can be seen 
here, indicating that the use of the low frequency catalogs is a good proxy for 
source detectability. 



9 Reliability 

The presence of an apparent source within the survey map does not guarantee 
that it is a real source that will be confirmed in a foUowup observation. The 
rate at which false positive detections are made by the source finding program 
should be taken into account. 

To estimate the reliability for each region of the sky within the AT20G 
a comparison was made between the sources that were found in the source 
detection process and a low frequency catalog that covered the same area of 
sky. Figure [T7] shows how the low frequency matching rate changes with the 
filtering metric. The rate at which sources are found within one of the three 
low frequency catalogs is used as the reliability and is denoted by TZ. At T 
values below -2.5 the fit approaches 0% and above -1.5 it reaches 100%. Neither 
of these two values are considered to be practical and the reported reliability, 
TZ, of a source is hmited to 5% <TZ < 98%. 

Once the relationship between and TZ had been established it was then 
possible to assign a reliability estimate to each source within the scanning 
survey catalog. Sources that were assigned a reliability oi TZ < 50% were 
removed from the final version of the scanning survey catalog. 
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Fig. 17 The fraction of sources that were matched with a low frequency catalog, as a 
function of the filtering metric. The red curve represents the relation between Ti. and T that 
is used to assign a reliability estimate to individual sources. The lower horizontal dotted 
line represents the spurious matching rate of 5%, whilst the upper solid line is the point at 
which half of the sources are matched. 

10 Survey Completeness 

In theory a comparison between a catalog and a list of all real objects gives 
an accurate measure of the catalog completeness, however an a 'priori list of 
all real objects does not exist. 

Due to the many calibration steps that were performed at the software level 
rather than the hardware level, and the relatively static performance of the 
correlator itself, it was possible to know a lot about how the system responded 
to a given source. Because of this it was possible to insert artificial sources 
into the data stream at the most basic level and then follow them through 
the calibration, map making, source finding and filtering process. Having this 
amount of control over the calibration and source detection process meant that 
it was possible to create an accurate measure of the scanning survey catalog 
completeness, and by extension, the follow-up survey catalog completeness. 

Survey completeness was estimated in a semi-analytical manner. In order 
to obtain a set of data that were consistent with the actual observations, whilst 
including as many of the detrimental factors as possible it was decided that 
the survey data would be used. Specifically, the regions between the detected 
sources were considered to be an excellent resource for the compilation of a 
data set that could be used to estimate the completeness and reliability of the 
survey. 

The injection of sources was achieved by inserting a scaled copy of the 
primary calibrator into the data stream as it was being processed. This method 
is sensitive to relative timing errors (ie. between scans) but not to absolute 
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timing errors. The absolute timing errors are detected via comparison of the 
final maps to known calibrators or the follow-up survey catalog positions and 
fluxes. The completeness estimate is thus assuming that the timing calibration 
is correct to first order. 

For each l'' x 10 - 15° map of the sky, a list of sources was produced 
that had a specific range of fluxes and a random distribution of positions. The 
source positions were constrained so that they were not within a primary beam 
width of an already detected source. At 20 GHz the density of sources detected 
is of the order 0.3 deg~^, and as the number of sources injected was around 500 
per map, the injected source positions were constrained so that the measured 
completeness would not be affected by excess amounts of confusion. In this way 
it was possible to create a map of the sky that contained all of the properties 
of the real observations that would otherwise be hard to characterise and 
reproduce in a simulated map. Such properties include but are not limited to: 
sensitivity variations, regions that were subject to patching (possibly multiple 
times), Tgys variations due to weather that were not captured in the calibration 
process, and data flagging due to solar system object^j. It should be noted that 
the completeness measurement is made using only the injected sources and 
thus deriving a source count model for the scanning survey catalog using the 
completeness and reliability measures can be done in a non-circular manner. 
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Fig. 18 The recovery rate for injected sources within the region —15° > <5 > —30°. Left: 
All sources. Right: Sources with TZ > 90%. The red bars are the fraction of sources that 
were recovered from the maps, with the blue line being a fit to equation [S] 

The survey completeness was determined by flrst measuring the rate at 
which sources were recovered from the injected maps as a function of flux. 
Figure flSl shows the recovery rate for all maps within the region —15° > S > 
—30°. To eliminate the effects of noise within the recovery rate measurements, 
the function in equation [3] was flt to the data, where S is the flux, a is the 
maximum completeness and b,c are shape parameters. 



^ The Moon, Sun and Jupiter were troublesome for bands above —30° declination 
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Once a fit to equation [3] had been obtained for eacli of the maps, the 
recovery ra te was converted to a source count using the source count model of 
Sadler et al . (2008) in equation[S]as shown in tabled The survey completeness 
at a flux S is then given by 



CiS) 



g{s)N{s) 

/f N{S) 



(4) 



where Smax is sufhciently large that the completeness is no longer changing. 
Although the source count model of ISadler et al. I (l2008l) has uncertainties of 
20 — 50% in the actual number of sources within a particular flux bin, the 
completeness estimate outlined above is sensitive only to errors in the relative 
number of sources within each flux bin. Thus only the errors in the slope 
of the source counts (a and /?) will affect the completeness estimates. The 
uncertainties in a and /3 are 3% or less, and when integrated over the range of 
50mJy to 500 mJy the total source count deviates by less than 9%, but this 
affects the calculation of equation 2] by less than 1%. 



parameter 



/3 



1.92 
2.78 
1.09 
28.9 



value 

Trm- 



N{S): 



2N, 



0.02 
+0.09 
0.07 
+0.13 
0.18 



Jy 



Table 2 The source count model at 20 GHz as reported by ISadler et al ] ll2008h . 



(5) 



Figure [TO] shows the corresponding completeness levels for the plots of 
figure UHl The effect of filtering is also demonstrated, and can be seen to 
lower the completeness of the catalog as the reliability is increased. Source 
that were within 1.5°of the Galactic plane are excluded from the follow-up 
survey catalog and are thus not included in the completeness calculations. 
The Galactic sources thus do not contaminate the completeness of the maps. 

The major contribution to the degradation of completeness seen in some 
maps in figure [20l is changes in sky coverage. The maps that have a com- 
pleteness below the average of their region have a different colour. In each of 
these cases the lower completeness is due to data that were flagged as poor but 
was unable to be recovered during patching observations or that was simply 
not observed due to weather and scheduling constraints. The scanning survey 
catalog is 98% complete at 70 mJy and 92% complete at 50 mJy. 
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Fig. 19 The level of completeness as a function of corrected survey flux density as measured 
by the injection and recovery of sources for the region —15° > S > —30°. The different 
coloured lines represent sources with a reliability estimate of 7?, > (0.05,0.5,0.75,0.9). The 
completeness is given by equation |4] 




Fig. 20 The completeness of the survey as measured via the injection and recovery of test 
sources, at 70mjy (Left) and 100 mJy (Right). 



11 Scanning Survey catalog 

The final version of tfie scanning survey catalog (SSC) is made of from sources 
with an SNR > 5 and a reliability of 50% or greater. Table [3] shows the first 
page of the scanning survey catalog, which contains a total of 7754 sources. 
The catalog format is as follows: 

- Columns 1&2: The source RA/DEC in J2000 co-ordinates 

— Columns 3&4: Galactic {l,b) 
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— Columns 5&6: Corrected survey flux at 20 GHz and the map rms noise 
(mJy) 

— Column 7: Reliability TZ 

— Column 8: Epoch of scanning observation. Refer to table [T] for the precise 
dates 

— C olumn 9: AT20G narn e, when the source is in the follow-up survey catalog 
of iMurnhv et~aD (l20inh 



The typical positional errors in a and 6 are ±15arc-sec. The flux scale has 
been corrected to be the same as the follow-up survey catalog. There may 
be variations in the individual fluxes due to source variability and the lower 
signal to noise of the survey maps. The AT20G names are determined by 
cross matching the follow-up survey catalog (FSC) and the scanning survey 
catalog (SSC). The nearest scanning survey catalog source within 120 arcsec 
is considered to be the correct identification. 



12 Summary 



The AT20G survey was made possible via the use of a wide band analog cor- 
relator in conjunction with an interleaved fast scanning strategy (or on-the-fly 
mapping). New programs were created that allowed for the offline calibration 
of the scanning survey data. The subsequent calibration scheme was accurate 
enough that the scanning data could be brought together to form convolved 
maps of the sky from which candidate source lists could be created. The offiine 
data processing techniques allowed artificial sources to be introduced into the 
data reduction pipeline which allowed for an accurate measurement of the com- 
pleteness of the scanning survey catalog, and b y extension, the com pleteness 



of the follow-up survey catalog as described in (|Murphv et all . l2010l ) 



The scanning survey catalog presented here has a lower flux limit and is 
more complete than the follow-up survey catalog, but achieves these goals at 
the cost of reliability. The follow-up survey catalog is ^^100% reliable, whilst 
the scanning survey catalog includes some sources with only 50% reliability. It 
is expected that the primary use of this c atalog will be t o provide foreground 
masks for CMB missio n s such as Planck (|Taubeii . 120051 ) and the South Pole 
Telescope (jStark et all . ll998V where completeness is more important than 
accurate fluxes or a high reliability. 
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Table 3 The first page of the scanning survey catalog. A pre-pubUcation version of the 
entire catalog is available at www.physics.usyd.edu.au/~haiicock/files/SurveyCatVl. 1 
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